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Abstract -We consider global quenches in the quantum XY chain in a transverse field and study 
the nonequilibrium relaxation of the magnetization and the correlation function as well as the 
entanglement entropy in finite systems. For quenches in the ordered phase, the exact results 
are well described by a semiclassical theory (SCT) in terms of ballistically moving quasiparticles. 
In the thermodynamic limit the SCT is exact for the entanglement entropy and its modified 
version following the method of Calabrese, Essler and Fagotti: arXiv: 1204-3911 is exact for the 
magnetization and the correlation function, too. The stationary correlation function is shown to 
be described by a generalized Gibbs ensemble. 



Introduction. — Recent progress of experimental 
work on ultracold atomic gases in optical lattices [THTO] 
allowed it to study the unitary time evolution of quan- 
tum systems in nonequilibrium situation, i.e. after a so- 
called global quench. This is usually achieved by a sudden 
change of the parameters of the quantum system within 
a time scale that is much shorter than the characteristic 
time the system needs to relax into a stationary state. 
The basic questions in this context are i) the dynamical 
characteristics of the relaxation process and ii) the pos- 
sible existence of a stationary state after long times. In 
three-dimensional systems fast relaxation into a thermal 
stationary state has been observed. In contrast to this, in 
quasi-one-dimcnsional systems the relaxation process has 
been found to be much slower and to lead to an unusual 
nonthermal stationary state [6]. This result has stimu- 
lated intensive theoretical work [111431] to clarify the ef- 
fect of integrability of the system on the relaxation pro- 
cess as well as on the nature of the stationary state. It 
is commonly expected that observables of nonintegrable 
systems effectively thermalize, which means that their sta- 
tionary state is described by a thermal Gibbs ensemble. 
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Numerical studies of different non-integrable systems are 
in favour of this expectation |12H21j . however some con- 
tradictory results indicate that the issue could be more 
complicated 22 24,30 . On the other hand in intcgrablc 
systems, due to the existence of integrals of motion, the 
stationary state is expected to be represented by a gen- 
eralized Gibbs ensemble (GGE) [T5], in which each mode 
corresponding to a conserved quantity is characterized by 
its own effective temperature. Results on integrable sys- 
tems are mainly collected on free fermion models, such 
as on the transverse Ising chain, for which several ana- 
lytical and numerical results have been recently obtained 
32 45]. Qualitative features of the relaxation process can 
be explained with a quasiparticle (QP) picture [T31H5] : 
The quench changes the total energy of the system by an 
extensive amount, which creates QPs homogeneously in 
space, moving ballistically with a constant velocity. Due 
to the conservation of momentum, the QPs are created in 
pairs with opposite velocity and these QP-pairs are quan- 
tum entangled. This QP-picture has been used to explain 
the time evolution of the entanglement entropy 35 , 47 and 
has been made quantitative with a semi-classical theory to 
predict also the relaxation of the local magnetization and 
correlation function [3111111113 ■ In intcgrablc systems like 
the transverse Ising chain QPs do not interact, giving rise 
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to a GGE description of the stationary state. In contrast 
to this, QPs in nonintegrable systems perform collisions 
and scattering, leading to thermalization. 

The scenario of thermalization or quasithermalization 
described above holds, strictly speaking, in the thermody- 
namic limit, when during relaxation of finite segments of 
the chain the rest of the system plays the role of the envi- 
ronment. However, the concrete experiments are usually 
performed on a few hundreds of atoms, therefore it is of 
importance to study finite size and boundary effects. Re- 
cent numerical investigation on the finite transverse Ising 
chain shows [MET] that i) the relaxation of the local mag- 
netization depends on the position and that ii) there is a 
quasiperiodic time dependence for long times. These fea- 
tures were explained and quantitatively well described by 
a semiclassical (SC) theory [UJ involving the reflection of 
the QPs at the open boundaries of the system. 

In this Letter we consider another magnetic model, the 
XY chain in a transverse field (denoted simply as XY 
chain), which is a generalization of the transverse Ising 
chain and which can also be mapped onto free fermions. 
This model has been studied in detail both in equilibrium 
[32149] and out of equilibrium [32135] after a global quench. 
In particular, the time evolution of the entanglement en- 
tropy has been obtained analytically in the large system 
limit [35]. Here we consider large finite chains and calcu- 
late the time dependence of the entanglement entropy, the 
local magnetization and the equal-time correlation func- 
tion by free fermion techniques [32][49]. These results are 
then compared with the prediction of our SC theory. We 
show that in the thermodynamic limit the SC theory pro- 
vides exact results for the entanglement entropy. For the 
magnetization and for the equal-time correlation function 
the SC theory can be modified along the lines of Ref. [IT] , 
so that it will become asymptotically exact [42][44], too. 

The remaining paper is organized as follows. First we in- 
troduce the model, present the elements of the free fermion 
solution and describe the basic ingredients of the SC the- 
ory. The main part contains our results, which are dis- 
cussed at the end. 

Model and the free fermion representation. 

The XY chain is defined by the Hamiltonian 
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in terms of the Pauli spin operators a x ' v ' z at site I = 
1 , . . . , L — 1 . Generally, we consider large finite chains of 
length L with open boundary conditions. The parameters 
< 7 < 1 and h > denote the strength of the anisotropy 
and the transverse field, respectively. The special case 
7=1 represents the transverse Ising model, and for h = 0, 
7 = the Hamiltonian reduces to the XX chain (see the 
equilibrium phase diagram in FigJTJ). 

We consider global quenches (at zero temperature), 
which suddenly change the parameters of the Hamilto- 
nian from 70, ho for t < to 7, h for t > 0. For t < the 



system is assumed to be in equilibrium, i.e. in the ground 
state of the Hamiltonian % with paramters 70 and ho, 
which is denoted by |$o)- After the quench, for t > 0, 
the state evolves coherently according to the new Hamil- 
tonian as |$o(*)) = exp(— i"Ht) |3>o)- Correspondingly the 
time evolution of an operator in the Hciscnberg picture is 
07 (t) = cxp (iHt) 01 cxp (—%Ht). 

We calculate the general two-point correlation function 

G xx (Mi; / a ,t a ) = ($ K(ti)af 2 (t 2 )\ $0) (2) 

in particular for t = t\ = t%, which is the equal-time cor- 
relation function denoted by C xx (h,fe)- For large sepa- 
rations it is C xx (h,h) = mh (t)mi 2 (t), where mi (t) is 
the local magnetization. In the initial state in the ther- 
modynamic limit (L — > 00) one has m; (0) > (mi (0) = 
0(1/L)) for ho < 1 (ho > 1) and at ho = 1 there is a 
quantum critical line, which belongs to the (transverse) 
Ising universality class [5D] for 7 > 0, i.e. mi (0) ~ L -1 / 8 . 

Using standard techniques [35], the Hamiltonian in 
Eq. (JJJ is expressed in terms of fermion creation and anni- 
hilation operators rp p and r\ v as 
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where the sum runs over L quasi-momenta, which are 
quasi-cquidistantly distributed in < p < ir, whereas for 
periodic chains we have \p\ < n. The energy of the modes 
is given by 
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and the Bogoliubov angle p diagonalizing the Hamilto- 
nian is given by tan(9 p = — 7 ship/ (h — coup). The cor- 
relation function in Eq.Q can be written as a Pfaffian, 
which is then evaluated through the determinant of an 
antisymmetric matrix. The local magnetization is calcu- 
lated in the form of the off-diagonal matrix element [501 
mi (t) = ($0 I of I $1)1 where |$i) is the first excited state 
of the initial Hamiltonian. In the free fermion description 
mi (t) it is given by a 21 x 21 determinant. 

The entanglement entropy between a block of I con- 
tiguous spins and the rest of the system is defined as 
Si = Tr,^ [ft log ,9;], where pi = Tr 4>i |$ )($ | is the re- 
duced density matrix, which evolves in time as pi(t) = 
exp(iHt)piexp(—iHt). The calculation of the dynamical 
entanglement entropy of free fermion models is described 
in Ref. [51]. 

The semiclassical theory. — As mentioned in the in- 
troduction, QPs are created after the quench at t = 0. For 
the XY chain these QPs are the free fermions described 
above. The wave packets formed by the free fermions move 
balUstically with a constant velocity ±t> p , which is ob- 
tained in the SC theory as 
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The position of the QP pairs at times t > can be easily 
calculated from their creation position and their velocity; 
in finite open chains the QPs are reflected at the bound- 
aries. 

We also need the creation probability /„ = 
fp (ho, 70; h, 7) of the QP pair. Here we make use of the 
fact that in a homogeneous system the QPs are created 
uniformly in space and that f p corresponds to the occu- 
pation probability of mode p in the initial state |$o), thus 
fp = ($o| VpVp l^o)- For the XY model it is expressed 
through the difference A p = p — P of the Bogoliubov 
angles as f p = \ (1 — cos A p ) with 

(cosp- ho) (cosp- h) +770 sin 2 p . . 

cos Ap = 7M^¥) • (6) 

We note that in open chains f p has a small position de- 
pendence near the boundaries, which to leading order can 
be neglected for long chains. 

To calculate the time evolution of the entanglement en- 
tropy, one observes that each pair of entangled QPs (and 
only those) with one partner moving within the block and 
simultaneously the other in the rest of the system, con- 
tributes an amount s p = — (1 — f p ) In (1 — f p ) — f p In f p to 
the entanglement entropy. Summing up the contributions 
of all QP pairs, one obtains the value of the entanglement 
entropy at the given time. 

To calculate the relaxation of the magnetization, one 
should notice that in the magnetic basis in Eq.([T]) the QPs 
represent kink-like excitations. As described first for ther- 
mal excitations in Ref. [51] and generalized afterwards for 
quantum quenches in Ref. [41] , a QP passing site I changes 
the sign of the local magnetization operator erf . If in a fi- 
nite system the same QP visits I several times, erf changes 
sign only if the number of visits is odd. Summing up 
the contributions of all QPs which have passed I before 
t, one obtains the local magnetization m; (i). Similarly 
the correlation function can be calculated. For a finite 
system the SC theory has been worked out for the trans- 
verse Ising model in Ref. [31]; for details we refer to this 
work. For quenches deep in the ferromagnetic phase an 
excellent agreement between the SC theory and the exact 
results is obtained. For quenches close to the critical point 
the agreement is less good due to the fact that in these 
cases the kinks are not sharply localized and the domain 
walls have a finite extent of the order of the equilibrium 
correlation length. In the thermodynamic limit this effect 
can be taken into account by using an effective occupation 
probability [3TJ 
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This result follows from an asymptotically exact evalua- 
tion of the correlation function by Calabrese, Essler and 
Fagotti for the transverse Ising chain |42[|33j . As a mat- 
ter of fact this evaluation works analogously for other free 
fcrmion models, such as for the XY chain. In the follow- 
ing we use this modified SC theory also for finite chains in 
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Fig. 1: (a) Equilibrium phase diagram of the XY model at T — 
0. (b) The arrows represent the six different quench protocols 
considered in this Letter. 



order to calculate the time evolution of the magnetization 
and the correlation function . 

Results. — We have performed quenches with six dif- 
ferent pairs of parameters (/io,7o) — > [h,l) as indicated 
in the phase diagram in Fig. [TJ Besides the relatively 
small quench in the ordered phase (I) we have performed 
larger quenches in the ordered phase (II and III) as well 
as quenches between the ordered and the disordered phase 
(IV and V) and a quench to the XX model (VI). 

Entanglement entropy. The dynamical entanglement 
entropy is calculated for a finite chain of length L = 
256. The block is represented by the first I spins with 
I = 32, 64, 96 and 128. For the six different quenches 
the results of the free fermion calculations arc shown in 
Fig. [2] together with the predictions of the SC theory. 
We observe an excellent agreement. In the SC calculation 
one should analyse the trajectory of the QP pairs as de- 
scribed below Eq.©. For free boundary conditions their 
contribution to the entanglement entropy is summarized 
as follows. If the QP pair is created at site j, then for 
t < T p = L/vp the pair contributes to the entropy within 
the interval t Pl i < t < t Pi i with = \l — j\ jv p and t Pl 2 = 
min [(I + j) /v p , (2L - (1 + j)) /v p }. For T p < t < 2T p this 
effective interval is t p ^ < t < t Pt 4 with t Pl 3 = 2T p — t p ^ 
and t P; 4 = 2Tp — t Pt \. For t > 2T p the process is repeated 
with a period of 2Tj,. 

As seen in Figd] the time evolution of the entropy starts 
linearly, reaching a ^-dependent saturation regime. After- 
wards there is a linear decrease and the process is repeated 
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Fig. 2: (Colour online) Time evolution of the entanglement 
entropy of the block of the first I sites of the chain with free 
boundary conditions after the six quench protocols in Fig[U 
Free fermion (SC theory) results are indicated by full (broken) 
lines (-1 = 32, - I = 64, -I = 96, -1 = 128). 

quasiperiodically as for the transverse Ising chain [53| with 
a period T pei iod = L/v max , where v max is the maximum 
group velocity max p {v p }. 

According to the SC theory formulated above, the en- 
tropy grows like 

i r 

Si(t)=t— J dpvps p , t<l/v max (8) 
in the L — > oo limit and its saturation value is given by 

i r 

S l{ t )= l 7T dps p , t^l/Vmax- (9) 

27T Jo 

These limiting values correspond to the exact results [35] . 

In passing we note that the SC calculation for periodic 
boundary conditions is straightforward and agrees even 
better wih the exact data. The reason for the better agree- 
ment is here the lack of reflections of the quasi-particles 
at the open boundaries. Since there arc no reflections the 
period is half the period for free boundary conditions. 

Local magnetization. We have calculated the time evo- 
lution of the local magnetization mi (t) for finite open 
chains of length L = 256 at the positions I = 32, 64, 96 
and 128. In Fig. [3] the exact free fermion results are 
compared with the predictions of the modified SC theory 
for the six different quench protocols indicated in Fig. [T] 



Fig. 3: (Colour online) Time evolution of the local magne- 
tization mi (t) after the six quench protocols in Fig[T] Free 
fermion (SC theory) results are indicated by full (broken) lines 
(- I = 32, - 1 = 64, -1 = 96,-1= 128). 

If the quench is performed between two paramter points 
within the ordered phase, the modified SC theory repre- 
sents an excellent description of the relaxation process, in 
particular for short times (t < T pcr i d/2). If the quench 
involves the disordered phase as well or the XX point the 
agreement is less good. However, the qualitative features 
of the time evolution of the magnetization is similar in 
all cases in Fig. [3] and it has the same form as found 
previously for the transverse Ising chain [321131]. The ex- 
ponentially decreasing initial part of mi (t) is followed by a 
quasi-stationary plateau and afterwards there is an expo- 
nentially increasing reconstruction part. Then the process 
is repeated quasiperiodically. 

In the L —> oo limit the local magnetization is given by 

mi (t) = mi (0) exp f J dpv p f p 8 (I - v p t) 

- l lJ o d P U e (V - o) (10) 

in the modified SC theory, which defines the quench- 
dependent relaxation time r mag and the correlation length 

£mag 

2 f n ~ 2 f 7 ' 

T mag = - / dpVpf p , Cmag = ~ / dp f p . (11) 

~ Jo ~ JO 

Correlation function. We have calculated 

the equal-time correlation function Cf x (r) = 



p-4 



Quantum relaxation in the XY spin chain. . . 





Fig. 4: (Colour online) Equal-time correlation function 
qxx ^ £ Qr g xec j ^j me i after the six quench protocols in Fig[T] 
as a function of the distance r. Free fermion (SC theory) re- 
sults are indicated by full (broken) lines (— t = 10, — t = 20, 
- t = 40, - t = 60, - 1 = 80, - t= 100). 



C xx (i^£±i,t ; )t j ; r = i,3 ; ...,i _ l between 

two sites symmetrically located between the boundaries 
of the system. For L = 256 and for t= 10, 20, 40, 60, 80 
and 100 the space-dependence is shown in Fig. @] for the 
same quenches as indicated in Fig. [TJ The comparison 
of the free fermion results with the predictions of the 
modified SC theory shows, also in this case, an excellent 
agreement if the quench is between two parameter points 
in the ordered phase. As for the magnetization, the 
agreement is less satisfactory if the quench involves also 
the disordered phase or the XX point. In the L — s- oo 
limit the correlation function is given by 



Cr(r) = Cg x {r)exp\-t- J dpv p f p 6 (r -2v p t) 

-l^J dp f p 6 (2v p t - r)\ , (12) 

in the modified SC theory, which is analogous to the for- 
mula for the local magnetization in Eq. (fl"0]) . Note that 
the time r corr and the length scale £ cor r of the correla- 
tion function is related to that of the magnetization as 
Tcorr = T mag /2 and £ CO rr = £mag as for the transverse Ising 
chain gj. 

Discussion. — We have studied the nonequilibrium 
relaxation dynamics of the quantum XY chain after global 



quenches. In particular, we have considered the time evo- 
lution of the entanglement entropy as well as the relax- 
ation of the magnetization and the equal-time correlation 
function in finite systems. The numerical results obtained 
with the free fermion techniques were compared with the 
predictions of the semiclassical theory. For the entangle- 
ment entropy an excellent agreement was found for all 
types of quenches and the SC results were shown to be 
exact in the L — > oo limit. For the magnetization and for 
the equal-time correlation function the SC theory provides 
a very good approximation if the quench is performed be- 
tween two points deep within the ordered phase. For two 
points in the ordered phase closer to the critical line we 
showed that one can modify the SC theory by introduc- 
ing an effective quasiparticle occupation number, see in 
Eq.©, which provides exact results in rmodynamic limit. 
In a finite system new phenomena such as reconstruction 
of the magnetization and the correlation function occur, 
which is explained in terms of reflected quasiparticles at 
the open boundaries of the chain or, for periodic boundary 
conditions, in terms of the recurrence of the quasiparticles. 

In the stationary state the length scale is given by 
£corr = £mag hi Eq. ([TT|) . which is to be compared with 
the equilibrium correlation length £t measured at a finite 
temperature T [32] 
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As for the transverse Ising chain we should define a mode- 
dependent effective temperature T e g (p) by equating the 
integrands 



f P 
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ln|l-2/ p | 
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2T cff (p) 



which leads to the relation 



f P = 



exp 



1 



(14) 



(15) 



Here on the right hand side there is the Fermi distribution 
function, thus the nonequilibrium occupation probability 
of the free fermion mode f p is equal to its equilibrium 
thermal occupation probability at T c ff (p) . We note that 
the same relation holds for the transverse Ising chain |41j 
and it is expected to be valid for free fermion models in 
general. We can thus conclude that after a global quench 
the quantum XY chain will not be described in terms of 
a thermal distribution, either, but reaches a stationary 
state in which the correlation function is described by a 
generalized Gibbs ensemble. 
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